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LONG-TERM  GOAL 

Derive  a  new  drag  law  and  roughness  length  relationship  for  the  coastal  zone. 

OBJECTIVES 

Our  primary  objective  is  to  isolate  the  influences  of  wave  age  and  fetch  on  the  drag  coefficient  and 
surface  roughness  length.  This  includes  examination  of  the  influence  of  internal  boundary  layer 
development  on  heat  and  momentum  fluxes  in  the  coastal  zone,  that  can  lead  to  large  deviations  from 
existing  similarity  theory.  The  second  main  objective  is  to  augment  the  wave  age  with  more  specific 
wave  properties.  The  final  objective  is  to  provide  the  data  sets  to  other  modeling  groups. 

APPROACH 

Our  initial  objectives  were  realized  by  first  quality-controlling  the  RASEX  data  and  intercomparing 
fluxes  between  different  levels.  Different  estimates  of  the  "observed"  roughness  length  using  the 
profile  and  eddy  correlation  methods  were  compared.  Using  the  observed  values  of  the  drag 
coefficient  and  roughness  height,  different  existing  relationships  were  tested  and  new  formulations  for 
the  transfer  coefficients  for  heat  and  momentum  were  developed.  The  analysis  has  been  extended  to  a 
much  larger  RASEX  data  set  outside  the  intensive  period,  which  includes  a  large  sample  of  offshore 
internal  boundary  layer  cases.  We  have  also  added  the  offshore  tower  data  from  the  University  of 
Uppsala. 

WORK  COMPLETED 

During  the  past  year,  we  have  evaluated  the  Charnock  relationship  and  simple  wave  age  models  using 
RASEX  and  the  Swedish  data  sets.  We  have  also  explored  the  relationship  between  the  aerodynamic 
roughness  length  and  atmospheric  stability. 

RESULTS 

The  Chamock  coefficient  relates  the  aerodynamic  roughness  length  to  the  surface  friction  velocity.  To 
our  surprise,  we  found  that  the  success  of  the  Charnock  relationship  is  almost  exclusively  due  to 
artificial  correlation  guarenteed  by  the  relationship  between  the  friction  velocity  and  the  roughness 
length  in  the  Monin-Obukhov  drag  law  (stability-modified  log  law).  In  fact  a  higher  order 
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parameterization  of  the  roughness  length  in  terms  of  the  friction  velocity  out  performs  the  Chamock 
relationship,  although  the  success  of  this  formulation  is  due  primarily  to  artificial  correlation  as  well. 
These  results  were  suggested  by  both  the  random  rearrangement  of  actual  values  of  the  input  variables 
(heat  flux,  wind  speed  and  momentum  flux)  and  by  specifying  a  gaussian  distribution  of  each 
independent  variable. 

We  have  found  that  the  aerodynamic  roughness  length  depends  substantially  on  atmospheric  stability. 
The  roughness  length  is  orders  of  magnitude  smaller  with  stable  conditions,  presumably  due  to 
supression  of  downward  momentum  transport  to  the  sea  surface  is  weaker.  This  stability  dependence 
also  leads  to  a  stability  dependence  of  the  Charnock  coefficient  (Figure  1).  For  these  analyses,  the 
RASEX  data  did  not  contain  a  sufficiently  wide  range  of  stability  to  assess  the  stability  dependence. 
We  therefore  supplimented  the  Swedish  data  sets  with  existing  data  sets  which  we  have  in-house. 


bulk  Richardson  number 


Figure  1.  The  dependence  of  the  Charnock  coefficient  on  the  bulk  Richardson  number  defined  in 
terms  of  the  temperature  and  wind  at  the  observational  level  and  the  sea-surface  temperature.  The 
data  sets  are  the  Swedish  8-m  data  from  an  offshore  tower  in  the  Baltic  ( red  solid),  LongEZ  aircraft 
data  from  approximately  15  m  during  SHOWEX  ( November  1997-black  dashed,  March  1999-blue 
solid,  November  1999-red  dashed)  and  TOAGA  COARE  NCAR  Electra  data  at  approximately  30  m 

(green  dashed), 
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All  the  data  sets  indicate  a  similar  decrease  of  the  Charnock  coefficient  with  increasing  stability  except 
for  the  TOAGA  COARE  data.  In  the  latter  data  set,  significant  stability  is  achieved  by  very  weak 
winds  rather  than  significant  air-sea  temperature  difference.  Furthermore,  the  number  of  cases  is  fairly 
small  because  of  the  minimum  wind  speed  criteria. 

However,  we  have  not  yet  reached  a  decision  on  whether  a  formulation  of  a  stability  dependent 
roughness  length  or  Chamock  coefficient  is  ready  for  use  in  models  where  no  information  on  wave 
state  is  available.  This  stability  dependence  is  different  between  weak  wind  and  moderate  wind 
conditions,  persumably  due  to  the  increased  influence  of  swell  with  weak  winds.  The  composited 
results  in  Figure  1  exclude  winds  less  than  4  m/s.  More  investigation  is  needed. 

We  have  also  developed  models  of  the  roughness  length  and  Charnock  coefficient  on  wave  age. 

Again,  the  resulting  relationship  is  largely  due  to  artificial  self-correlation.  We  are  currently  studying 
the  influence  of  swell  propagation  with  respect  to  the  mean  wind. 

IMPACT/APPLICATION 

The  previous  success  of  the  Charnock  formulation  in  the  literature  may  have  also  been  primarily  due  to 
artificial  self-correlation.  This  does  not  suggest  that  common  use  of  the  Chamock  formulation  in 
numerical  models  should  be  terminated  but  rather  that  existing  physical  interpretation  of  the  Charnock 
formulation  cannot  be  justified.  In  addition,  the  common  use  of  the  Charnock  relationship  to  reduce 
the  drag  coefficient  to  neutral  values  is  not  justified  and  may  lead  to  misleading  results. 

RELATED  PROJECTS 

Work  on  an  ONR  grant  entitled  "Spatial  Variations  of  the  Wave,  Stress  and  Wind  Fields  in  the 
Shoaling  Zone"  (N00014-97-1-0279)  has  completed  the  observational  work  and  is  concentrating  on 
data  analysis.  This  program  studies  spatial  variations  in  the  coastal  zone  using  the  FongEZ  research 
aircraft  and  ground  based  sonic  anemometers  at  the  shore  and  on  the  pier. 

SUMMARY 

We  have  analyzed  a  variety  of  offshore  data  sets  to  study  of  fluxes  of  heat  and  momentum  between  the 
sea  surface  and  the  atmosphere  in  the  coastal  zone.  We  have  found  that  the  vertical  structure  of 
atmospheric  offshore  flow  is  much  more  complex  than  previously  thought,  especially  in  the  case  of 
advection  of  warm  air  from  land  over  cooler  water.  We  have  developed  new  formulations  for 
modeling  sea  surface  fluxes.  Modification  of  computer  models  to  predict  the  atmospheric  structure  in 
the  coastal  zone  requires  more  work  for  cases  of  complex  seas  and  warm  air  over  cooler  water. 
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